Damping of phase fluctuations in superfluid Bose gases 
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Using Popov's hydrodynamic approach we derive an effective Euclidean action for the long- 
wavelength phase fluctuations of superfluid Bose gases in D dimensions. We then use this ac- 
tion to calculate the damping of phase fluctuations at zero temperature as a function of D. For 
D > 1 and wavevectors |fe| <C 2mc (where m is the mass of the bosons and c is the sound ve- 
locity) we find that the damping in units of the phonon energy Ek = e[As| is to leading order 
Jk/Ek — An(kQ /2-Kp)(\k\/ko) 2D ~ 2 , where p is the boson density and ko = 2mc is the inverse healing 
length. For D — > 1 the numerical coefficient Ad vanishes and the damping is proportional to an addi- 
tional power of j ft | /ko ; a self-consistent calculation yields in this case y k /E k = 1.32 (k /2np) 1/2 \k\/k . 
In one dimension, we also calculate the entire spectral function of phase fluctuations. 
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I. INTRODUCTION 

It is well knowrP^ that the perturbative treatment of 
fluctuation corrections to Bogoliubov's mean-field the- 
ory^ f or interacting bosons is plagued by infrared diver- 
gencies, which appear at zero temperature for dimen- 
sions D < 3, and at finite temperature for D < 4. The 
physical origin of these divergences is the coupling be- 
tween transverse and longitudinal fluctuations^. As a 
consequence, the anomalous part of the single-particle 
self-energy 2^(0) at vanishing momentum and frequency 
is exactly zerc0, whereas Bogoliubov's mean-field theory 
predicts that Sa(0) is finite. To recover the exact result 
^U(O) = diagrammatically, infinite orders have to be 
re-summed using non-perturbative methods, such as the 
renormalization groupPEGJ 

If one is interested in long- wavelength and low-energy 
properties of t he sy stem, Popov's quantum hydrody- 
namic approactpJ^^ offers an alternative parametrization 
of the fluctuations which does not lead to infrared diver- 
gencies. In this approach one separates the low-energy 
from the high-energy modes and treats the low-energy 
sector within a gradient expansion for the phase and am- 
plitude fluctuations. This hydrodynamic approach can 
also be used to study interacting bosons in one spa- 
tial dimension, where strong fluctuations prohibit the 
formation of a Bose-Einstein condensat e ^ 12 * 13 ! , although 
the groundstate is superfluid. In fact, in one dimension 
the weak coupling expansion of thermodynamic quanti- 
ties obtained within the hydrodynamic approach agrees 
with exact results for the Lieb-Liniger modef^ up to 
the second order in the relevant dimensionless interaction 
parametei^. On the other hand, the single-particle spec- 
tral function and the dynamic structure factor (spectral 
function for density fluctuations) of interacting bosons in 
one dimension ha ve rec ently been shown to exhibit alge- 
braic singularities! 15 * 16 ^ In principle it should be possible 



to reproduce these singularities within the hydrodynamic 
approach, but this requires a non-perturbative treatment 
of the interactions between amplitude and phase fluctu- 
ations which is beyond the scope of this work. 

Here we shall use the hydrodynamic approach to cal- 
culate the damping of phase fluctuations in low dimen- 
sional Bose gases. In one dimension we also calculate the 
entire spectral function of phase fluctuations and show 
that in the vicinity of the phonon peaks it has approxi- 
mately Lorentzian line-shape, with on-shell damping pro- 
portional to k 2 for small wavevectors k. We also elabo- 
rate on the relation between the fc 2 -scaling of the damp- 
ing in D = 1 and the Beliaev damping of the phonon 
mode in superfluid Bose gases, which in D > 1 is known 
to scale as |fc| 2£>_1 for small wavevectors^. 

II. EFFECTIVE ACTION FOR PHASE 
FLUCTUATIONS 

According to Popo\fl^the long- wavelength asymptotics 
of correlation functions of interacting bosons can be ob- 
tained from an effective long-wavelength hydrodynamic 
action involving a phase field tp(r,r) and a conjugate 
density field p(r,r). These are slowly varying functions 
of space r and the imaginary time r, and are defined by 
writing the slowly varying part of the fundamental boson 
field as 

yj{r,T) = y^7)e 1 ^^. (1) 
Setting p(r, t) = p a + <r(r, r), where 

Pa = j d D r J drp(r,T) (2) 

is the spatial and temporal average of the density field, 
and expanding the effective action of the slow part of the 



2 



boson field to second order in the gradients, we obtain the 
hydrodynamic Euclidean action for the slowly varying 
phase and amplitude fluctuation;^ 



S[y,<?] = -/3Vp(p,p ) + S 2 [cp,a], 



(3) 



where /3 is the inverse temperature, V is the volume of 
the system, and p(p,po) is the pressure as a function 
of the chemical potential p and the average density po, 
and S 2 [<p, c] contains fluctuation corrections up to second 
order in the derivatives, 



S 2 [tp,<r] 







dr I d D r 



-iPu. Po <rd T cp-p PoPo - 



2m 
2 , (Vcr) : 



8mpo 



2m 



(4) 



Here m is the mass of the bosons and the coefficients 
Phi Phhi Ppopo an d Pp Po are the partial derivatives of the 
pressure p(p, po) of a homogeneous system with chemical 
potential p and average density pq. The last two terms 
on the right-hand side of Eq. Q represent the kinetic 
energy of the slowly oscillating part of the boson field. A 
simple approximation for the pressure is^ 



P(M)Po) ~ PPo 



u 2 



Uq 



[(A) - P? 



(5) 



where Uq is the two-body interaction at vanishing exter- 
nal momenta, and p = p/uq is the value of the fluctu- 
ating variable po at the saddle point of the functional 
integral. In the thermodynamic limit and at zero tem- 
perature we may identify p with the total density of the 
bosons. Eq. §5§ implies the following estimate for the 
relevant partial derivatives of the pressure, 



Pp~ P = 
Pftfi ~ 0, 
Pp.po ~ 1> 
PpoPo ~ ~ u 



p/u , 



(6a) 
(6b) 
(6c) 
(6d) 



The above hydrodynamic action describes long- 
wavelength fluctuations at length scales larger than 
some cutoff scale 1/A . In momentum space we should 
therefore impose an ultraviolet cutoff Ao on all integra- 
tions. In the weak coupling regime a reasonable choice 
of the cutoff is the inverse healing length Ao = 2mc, 
where c is the sound velocity defined below. 

Introducing the Fourier transform of the fields in 
momentum-frequency space, 



<p(r, T )= I e^™W 



K 



(7a) 

a(r,r)= I e^'™W, (7b) 



K 



that the gradient contribution Q to the hydrodynamic 



action can be written as 



-(- 



K X JK 2 J K 2 



PpoPo 

Sk 1 +k 2 +k 3 ,o 



4m po 
ki ■ k 2 



J(T-K&K 



(8) 



Here K — (k, iuS] is a collective label for momenta k and 
bosonic Matsubara frequencies iuj, the integration sym- 
bols represent J K = {fiV)^ 1 J2k 2~L' an d the normaliza- 
tion of the delta-symbols is 8k,K' — PVSk_k'Suj,uj' where 
the (5-symbols on the right-hand side are Kronecker- 
deltas. Since the hydrodynamic action ([8| is quadratic 
in the amplitude field a, we may carry out the functional 
integration over the cr-field, 



= / V\<j\e 



-Sa[<p,tr] 



The effective action of the phase field is 
1 



(9) 



Se{[[tp] 



G {K)(p_ K (p K 



1 

3! 



A" 



K ± JK 2 JK 3 



^>K 1 +K 2 +K 3 ,0 



xI ( K 1: K 2, K 3 )<p Kl tp K2 ip K3 

4} I I I I 3k 1 +k 2 +k 3 +k 4 ,o 



K1 JK 2 JK, JK. 



x T 4) (Ki ,K 2 ,K 3 , K 4 )ip Kl ip K2 ip K3 ip Ki , (10) 
where the inverse Gaussian propagator of the phase field 



Ppp 



P 2 



k 2 

Po T impa 



-Ppo 



(11) 



and the properly symmetrized three-point and four-point 
vertices are 



r,(3W ry- ts- s k x k 2 
T '(Ki,K 2 ,K 3 ) = -p WQ 



Ppopo + impo 

+(K 2 ^K 3 ) + (K 1 ^K 3 ), (12) 
(fei ■ k 2 )(k 3 ■ k 4 ) 



Oi+fe 2 ) 



m y Ppopo + 4mp 

+{K 2 ^ K 3 ) + {K 2 ^ K 4 ). (13) 
Note that the non-Gaussian contributions to the effective 



hydrodynamic action (10) of the phase fluctuations are 
generated by the term (Vy) 2 cr/(2m) associated with the 
coupling between amplitude and phase fluctuations in our 
original hydrodynamic action Q. 



III. DAMPING OF PHASE FLUCTUATIONS IN 
DIMENSIONS D > 1 

Within the Gaussian approximation we obtain the en- 
ergy dispersion E k of the phase fluctuations from the 
condition Gq E k + i-q) — 0. Approximating the pres- 
sure derivatives by Eqs. (|6aMdl we obtain 



G (K) 



u (l + k 2 /k 2 ) 



where E k is the Bogoliubov dispersion, 



E k =c\k\^/l + k 2 /kl 
Here the sound velocity is given by 



m 



and the inverse healing length 



(14) 



(15) 



(16) 



(17) 



marks the crossover from the linear regime of a sound-like 
dispersion to the quadratic regime of quasi-free bosons. 
Note that the bare coupling can be written as 



mc 
P 



(18) 



which in one dimension has units of velocity. In fact, 
in D = 1 the dimensionless ratio uq/c = mc/ p can be 
identified with the usual Lieb-Liniger parameter^ which 
is the relevant dimensionless interaction strength. 



The interactions in our effective action ( 10 1 give rise 



to a momentum- and frequency dependent self-energy 
T,(K), so that the true inverse propagator of the phase 
fluctuations is 



G- 1 (K) = Gv 1 (K) + X(K] 



(19) 



The renormalized energy dispersion Ef, of the phase 
mode and its damping "f k are given by the solutions of 
G~ 1 (k,E k + vy k ) = 0. To first order in the quartic ver- 



tex Fg 4 -' and to second order in the cubic vertex T^' the 
self-energy is T,(K) = T, r (K) + T, 2 (K), where 

Ei(A-) = \ J G (K')T^(K',-K',K,-K), 

MK) = ~I G (K')G (K' + K) 
1 Jk' 



(3) 



(20) 



xT^iK, -K - K',K')T^'(-K', K + K' , -K). (21) 

The corresponding Feynman diagrams are shown in 
Fig. [T] To lowest order in perturbation theory the damp- 
ing of the phase mode is given by 



.(3) 



Ik 



u (l + k 2 /k 2 ) 
2E k 



ImE 2 (k,E k + ir)), (22) 




1 

'2 



FIG. 1. These Feynman diagrams represent the first two per- 
turbative corrections to the self-energy of the phase fluctua- 
tions, see Eqs. (20] 21 1. The soli d lin es represent the Gaussian 
propagator Go(K) given in Eq. (141, while the black triangles 
and the black square denote the symmetrized three-point and 
four-point vertices defined in Eqs. \12\ |13|). 



where rj > is infinitesimal. Substituting Y.2(K) from 

Eq. J2lj) into Eq. (1221) and using Eqs. (fill) and Jl2l for 
1 — r . 1 — r 1 — 1 1 — r 

Go and T , we obtain for |fc| <C k after straightforward 
algebra 



7fc 



with 



71"»0 

16m 2 



^-^S(E k - E k , - E k _ k ,)W k , k >, (23) 

(ZTT) 

[k 2 - k' 2 } 2 [k 2 - (k - k') 2 } 2 



E k - k > E k < 
[k' 2 - (fe - k') 2 ] 2 

E k 



(24) 



Taking into account that the function W k k i is multiplied 
by 5{E k — E k i — E k - k i) 1 we may substitute under the 
integral sign for small momenta 



ir. 



9 



k.k' 



k\\k'\\k - k'\ 



(25) 



The fe'-integration can now be performed using D- 
dimensional spherical coordinates. For small external 
momentum k the loop momentum k' is almost parallel 
to k so that we may approximate^ 



5{E k — E k i — E k . 



|fc|-|fc'l 



V^clfellfe' 



(26) 



where i9 is the angle between k and k! . We finally obtain 
in D dimensions 



^=A D k l 
E k 



27rp y fcp 



Ifcl 



2(D-1) 



(27) 



where the numerical coefficient Ajj can be expressed in 
terms of T-functions as follows, 



A D = 



T(D) 



23£>r(^=i)r(£> • 



(28) 



4 



0.15 



0.1 - 



0.05 



intermediate states in our perturbative self-energy (21 1, 



we replace the Gaussian propagators on the right-hand 
side by the exact propagators G{K) of the phase mode, 
for which we use the spectral representation 

duo' B(k,u') 



2 
D 



FIG. 2. Graph of the numerical coefficient Ad defined in 



G(fc,tw)= / (29) 

J_ OD 2?r w u + uj 2 

where the spectral function 

B{k,uj) = 2ujlmG(k,uj + ir]) (30) 

is real and positive. Retaining only the imaginary part of 
the self-energy, we find after analytic continuation that 
for frequencies close to ±-Eft the spectral function can be 
approximated by 



Eq. (281 as a function of D for 1 < D < 3. 



B(k,w) 



u j(k,uj) 



(\co\-E k ) 2 + 1 2 (k^y 



(31) 



A graph of Ad as a function of the dimensionality D of 
the system is shown in Fig. [2j In D = 3 and D = 2 we 
obtain A 3 = jjg w 0.0187 and A 2 = *f « 0.108, while 
A n ~ tt(D — 1) — > for D -» 1. In three dimensions 



where the damping function 7(fc,w) satisfies the integral 
equation 

sgnw f dfc' r\"\ 
2^ _ . 

fc' 2 ! 2 



Eq. ( p7| agrees with the well-known Beliaev damping of 
the phonon mode in a Bose condensate^. Beliaev damp- 
ing in D — 3 and D — 2 has recently been re-derived in 
RefP^l using Popov's hydrodynamic approach; however, 
these authors did not integrate out the amplitude fluctu- 
ations, which renders the algebra more complicated than 
in our approach based on the effective action of phase 
fluctuations. The fact that for arbitrary D > 1 Beliaev 
damping scales as |fe| 2 ^~Mias been pointed out previ- 
ously by several authors^ 17 * 18 !. 



xB(k 
[fc 2 - 



16m 2 uo 
- k'.co - 



dLj'B(k',u>') 



'0 

ffc 2 - 



(k 



\UJ\ - UJ' 

[fc' 2 - (fc - fc') 2 ] 2 



(32) 



To solve this non-linear integral equation, we make the 
ansatJ^Sl 



y(k,u) = 7ft/ 



Ik 



(33) 



IV. 



PHASE FLUCTUATIONS IN ONE 
DIMENSION 



where the on-shcll damping is assumed to be of the form 
7fe = /ol^l"; with some exponent a. The dimensionless 
function f(z) is normalized such that /(0) = 1, so that 
the dimensionful constant /o determines the strength of 
the on-shell damping. The function f(z) is expected to 
be strongly peaked to z = and to decay as a power 
law for \z\ ^> 1. After substituting the ansatz (33) into 



uctuations in one dimension, because 
the coefficient Ad vanishes for D — >• 1. The problem is 



Obviously, Eq. (27} cannot be used to estimate the Eq . we may scale out the fc-dependence by intro- 
damping of phase fl 



that in the derivation of Eq. (27 1 we have inserted bare 



Green functions in the loop integration, which in D = 1 is 
not accurate enough to obtain the damping of the phase 
fluctuations. A similar problem arises in the calculation 
of the damping of the excitations of a clean Luttinger 
liquid, which has been studied by SamokhirPS by means 
of a self-consistent perturbative calculation taking the 
damping of intermediate states into account. Although in 
this case the spectral function is known to exhibit a non- 
Lorentzian line-shape with algebraic singularities,^ the 
overall width of the spectral function can be estimated 
correctly with this method. 

Let us now use the method proposed by SamokhirPSl 
to calculate the damping of phase fluctuations in the 
one-dimensional Bose gas. In fact, we shall go beyond 
Samokhin's work and calculate the entire spectral line- 
shape of phase fluctuations. To include the damping of 



ducing dimensionless integration variables x — fe'/fc and 
y = (uj' — Ey ) /7ft' . It is then easy to see that our ansatz 
is only consistent if a — 2. The constant /o is then given 
by 



/o = 



3/M/I k Q 



4m V 27r j o' 

where the function f(z) satisfies the integral equation 



(34) 



/(*) = 



(35) 



with the non- linear functional I z [/] given by 



Jo 1 ~ x J_, 



dy 



f(y) 



y 2 + P(y) 



f ( (z-y)x 2 \ 
J { (1-s) 2 J 



I (l-x) 2 J + J V (l-x) 2 I 



(36) 
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-1 1 

(M - E k )/ lk 



FIG. 3. The solid line is the numerical solution of the integral 
equ ation (351. The dashed line represents the approximation 
(37 1 with J2 = 0.22, which is reasonably accurate for \z\ < 5. 



The dotted line represents the interpolation ( 38 1 with /i = 



0.58. The logarithmic plot in the inset shows that for large 
| z | the solution of the integral equation (35 I vanishes as l/|z|. 



The dashed-dotted line in the inset is the curve 1 / (0.58 j 2; j 



FIG. 4. Normalized spectral function "/kB(k,oj)/uo of phase 
fluctuations of the one-dimensional Bose gas as a function 
of frequency for fixed k = 0.5 fco and Lieb-Liniger parameter 
uo/c — 0.1. The on-shell damping is in this case 7fc/(cfco) = 
0.059. The dashed line represents a fit to a Lorentzian with 
on-shell damping 7fc. 



The integral equation ( 35 ) can easily be solved numeri- 



cally. In practice we obtain convergence for any reason- 
able initial guess for the function f(z). It turns out that 
the ansatz 



/(*) 



1 



1 + fiz 2 



(37) 



with f 2 = 0(1), does lead to a rather fast convergence 
after a few iterations. The solution of the integral equa- 
tion ( 35 1 is represented by the solid line in Fig. [3] In 
fact, for /2 = 0.22 the ansatz (37 1 is already a reasonable 



approximation to the solution of Eq. (351 in the regime 
\z\ < 5. Because the quadratic z-dependence of the ex- 



act solution for small z is correctly described by Eq. ( 37 1 , 



this ansatz describes the spectral function in the vicin- 
ity of the quasi-particle peaks quite accurately. On the 
other hand, as shown in the inset of Fig. [3j for large \z\ 
the numerical solution of Eq. (35 1 decays as l/|z|, which 

The tails 



is not correctly described by our ansatz ( [37] ) 
of the spectral function are therefore better described by 
the interpolation formula 



/(*) 



1 



(38) 



which for /1 ps 0.58 has the correct asymptotics for large 
\z\, but is less accurate than Eq. ( |37[ ) for small z. 

Given our numerical solution f(z) of the integral equa- 
tion ( 35 ) , we obtain 



7o[/]«0.78, 



and hence 




(39) 



(40) 



We conclude that for small wavevectors the on-shell 
damping of the phase mode in the one-dimensional Bose 
gas in units of its energy Ej~ ps c\k\ can be written as 



lk_ 



1.324 



2np k 



(41) 



Note that the dimensionless ratio k /2np = mc/-K p = 
uq/ttc can be identified with the Lieb-Liniger parameter 
divided by 7r. Keeping in mind that in the derivation 
of Eq. (41) we have neglected vertex corrections, we ex- 



pect that the prefactor in Eq. (41) is accurate as long 
as the Lieb-Liniger parameter uq/c is small. Comparing 
Eq. (41 ) with the corresponding expression (27) in D > 1 



we see that in one dimension the damping involves an 
addition factor of |/c|/fc ; however, in D = 1 the prefac- 
tor is proportional to the square root of the Lieb-Liniger 
parameter uq /c, whereas in D > 1 it is linear in the cor- 
responding dimensionless parameter fcf /27rp. 



Since the solution of the integral equation ( 35 ) gives 



the entire scaling function f(z) in Eq. (33), it is now 



easy to obtain the momentum- and frequency depen- 
dent spectral function B{k,u) of phase fluctuations of 
the one-dimensional Bose gas. The result is plotted in 
Fig. [4] Obviously, for frequencies not too far away from 



the central peak 



Ek\ < 27fc) the line-shape can be 



approximated by a Lorentzian, but outside this regime 
the spectral function decays faster. Using the fact that 
f(z) ~ (AM) -1 for large z we find that the tails of the 
spectral function are 



u j(k,uj) 



1 



E k y fx \\lu\ - E k \ 3 ' 



(42) 



which decays faster than a Lorentzian by a factor of 
lk/\H-E k \. 



() 



V. SUMMARY AND CONCLUSIONS 



In summary, we have derived an effective action de- 
scribing the dynamics of low-energy and long- wavelength 
phase fluctuations of superfluid bosons. Using this ac- 
tion, we have then calculated the leading momentum de- 
pendence of the damping of the phase fluctuations in 
arbitrary dimensions. For D > 1 a simple perturba- 
tive calculation yields the usual Beliaev damping, which 
scales as |fe| 2D_1 in D dimensions. For 13 — » 1 the pref- 
actor of |fe| 2£>_1 vanishes, and the damping is propor- 
tional to k 2 . We have obtained this result by taking the 
damping of the intermediate states in the loop integra- 
tion self-consistently into account. In one dimension, we 
have also calculated the spectral function of phase fluc- 
tuations, which has a Lorentzian line-shape for frequen- 
cies close to the quasi-particle peaks associated with the 
sound mode, but for larger deviations from the peaks 
decays faster than a Lorentzian. 

Since the vertices of the effective action for the phase 
fluctuations vanish for zero wavevectors or frequencies, 
we believe that higher orders in perturbation theory do 



not qualitatively modify our results. In particular, in 
one dimension the spectral function of phase fluctuations 
does not contain any algebraic singularity, in contr ast to 
the spectral function of the amplitude fluctuation d 15 * 21 l 
We are not aware of experimental methods to directly 
measure the spectral function of phase fluctuations, so 
that we cannot compare our result for the spectral line- 
shape with experiments. However, for non-perturbative 
calculations of the single-particle Greens function of su- 
perfluid bosons one usually assumes that the Gaussian 
approximation is sufficient to calculate the propagator of 
the phase fluctuations^^. Our results imply that the 
Gaussian approximation is indeed well justified in this 
case, because the damping of the phase mode is small, 
so that in the superfluid state the phase fluctuations can 
propagate as well-defined quasi-particles, even in one di- 
mension. 
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